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Abstract. In second order perturbation theory different definitions are known of gauge in- 
variant perturbations in single field inflationary models. Consequently the corresponding gauge 
invariant cubic actions do not have the same form. Here we show that the cubic action for one 
choice of gauge invariant variables is unique in the following sense: the action for any other, 
non-linearly related variable can be brought to the same bulk action, plus additional boundary 
terms. These boundary terms correspond to the choice of hypersurface and generate extra, 
disconnected contributions to the bispectrum. We also discuss uniqueness of the action with 
respect to conformal frames. When expressed in terms of the gauge invariant curvature per- 
turbation on uniform field hypersurfaces the action for cosmological perturbations has a unique 
form, independent of the original Einstein or Jordan frame. Crucial is that the gauge invari- 
ant comoving curvature perturbation is frame independent, which makes it extremely helpful 
in showing the quantum equivalence of the two frames, and therefore in calculating quantum 
effects in nonminimally coupled theories such as Higss inflation. 
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1 Introduction 

The detection of non-Gaussianities in the Cosmic Microwave Background radiation (CMB) would 
provide a wealth of information. From the point of view of inflationary physics, it would con- 
strain the correlation functions of primordial fluctuations. Since these correlation functions can 
be explicitly found for one of the many inflationary theories, non-Gaussianity provides a pow- 
erful tool to constrain the parameters in these theories and discriminate between models in the 
inflationary zoo. 

One of the simplest inflationary models is that of a single scalar field in a slow-roll scenario. 
In his seminal work [1] Maldacena found the third order action for inflationary perturbations 
and showed that non-Gaussianities are too small to be observed. This has been extended to 
more general scalar theories in [2] and [3]. 

Importantly, the cubic action and 3-point functions have been derived for £, the curvature 
perturbation. It is well-known that £ itself is a gauge dependent quantity, in the sense that it 
is not invariant under reparametrizations of coordinates. Only a gauge invariant perturbation 
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can be called physical. On the other hand one can work with gauge fixed quantities, although 
one must be careful with respect to gauge artifacts. 

In this work we set out to find the cubic action for gauge invariant cosmological pertur- 
bations. We show how to do this for different gauge invariant variables at second order and 
that the gauge invariant actions reduce to the ones in [1] in the spatially flat or uniform field 
gauge. In general, the gauge invariant cubic actions for different variables are different, but 
we show that the bulk part of the gauge invariant action coincides for different variables. In 
that sense the evolution of non-Gaussianity is unique. The difference between the actions lies in 
boundary terms, which are associated with the choice of hypersurface. They generate additional, 
disconnected parts of the bispectrum. 

In the second part of this work we discuss uniqueness from the point of view of conformally 
related frames. As is well known, nonminimally coupled actions, or Jordan frame actions, can be 
written into a minimally coupled form, the Einstein frame, by field dependent redefinitions of the 
metric and scalar field. Thus in principle nontrivial Jordan frame results can be obtained from 
well-known Einstein frame results by redefining fields, which makes the frame transformation 
a very powerful tool. The situation becomes complicated at the level of perturbations. The 
problem is that the perturbations of the metric and scalar field in the Einstein frame are not 
equal to those in the Jordan frame, precisely due to the field dependent field redefinition. The 
situation is very similar to the gauge problem for perturbations. 

In this work we point out that it is possible to construct frame independent cosmological 
perturbations, which are very useful to relate Einstein frame results to Jordan frame results. 
Here is where the second aspect of uniqueness comes in: the cubic action for cosmological 
perturbations takes a unique form, independent of the frame, provided one makes use of precisely 
that variable that is the same in either frame. As it turns out, this variable coincides with the 
gauge invariant comoving curvature perturbation. 

A short outlook: in Sec. 2 we define the single scalar field action and its perturbations, 
in Sec. 3 we construct different gauge invariant variables at second order and in Sec. 4 we 
construct the gauge invariant action for cosmological perturbations and discuss uniqueness. In 
the second part, Sec. 6 we discuss different frames and perturbations in those frames, and show 
in what aspect the action is unique. 

2 Action and perturbations 

The action under consideration is the Einstein-Hilbert action for General Relativity plus a 
minimally coupled scalar field 



Here we use natural units in which 8irG = 1, c = 1 = K. This so-called Einstein frame action is 
manifestly covariant, that is, it is invariant under spacetime coordinate reparametrizations. It is 
possible to study the background field equations for (2.1). Taking as the background spacetime 
the homogeneous and isotropic FLRW metric leads for example to the Friedmann equations in 
presence of a time dependent scalar field, which give rise to inflationary solutions in certain 
regimes. 

It is also possible to study perturbations of the metric and the scalar field in the action (2.1). 
A perturbation of a quantity is defined as the difference between the quantity in the physical 
spacetime and the quantity in the background spacetime. In order to compare these quantities 
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one has to choose a mapping between the physical and background spacetimes. The gauge 
freedom in General Relativity is the freedom in choosing a mapping. As a consequence the 
perturbations themselves depend on the choice of mapping. They are in other words gauge 
dependent. 

Gauge dependence is in principle problematic, since physical results should not depend on 
reparametrization of coordinates. However, we know for a fact that the original, unperturbed 
action (2.1) is explicitly covariant. Therefore, it may be possible to write the perturbed action 
in a manifestly covariant way. Indeed, this can be achieved using gauge invariant cosmolog- 
ical perturbations [4]. Gauge invariance here is synonymous to covariance or diffeomorphism 
invar iance. 

A convenient method to deal with the gauge dependence in the action is to use the ADM 
formalism [5] with line element 

ds 2 = -N 2 dt 2 + g ij (dx i + N i dt){dx j + N j dt) . (2.2) 

Geometrically, spacetime has been sliced up in spatial hypersurfaces whose geometry is described 
by the spatial metric gij . The slicing and threading of spacetime is described by the lapse function 
N and shift functions N{. The action (2.1) with the ADM metric (2.2) becomes 

S = - J d 3 xdtN^g-^R^ + N- 2 (E ij ' - E 2 ) +N~ 2 (d <S> - A^) 2 - g^d^dj® - 2V($) j , 

(2.3) 

where 

Eij = \ (dogij - ViNj - VjNi) 
E = g^Eij , (2.4) 

and i?( 3 ) is the spatial scalar curvature computed from (spatial derivatives of) gij alone. From 
Eq. (2.3) it is clear that N and Ni are non-dynamical fields, and moreover they are gauge 
dependent. In fact, in a Hamiltonian formulation they appear as Lagrange multipliers. Solving 
for these fields in the action (2.3) removes the unphysical degrees of freedom. In the end the only 
dynamical perturbations out of the 7 degrees of freedom in g^ and are one scalar degree of 
freedom and a graviton. The remaining 4 degrees of freedom are nondynamical and are actually 
the solutions of the constraint equations [6] . The dynamical degrees of freedom are indeed gauge 
invariant. 

The perturbed action is calculated by inserting 

$ = (j)(t) + ip 

9ij = a{t) 2 (^e 2C + ^dtdfh + ^%^) + h ¥^j 
N = N(t) (l + n) 

Ni = a(t)N(t)(-di,s + nj) , (2.5) 
a 

where £ is the curvature perturbation and all perturbations carry a temporal and spatial depen- 
dence. We have made use of the scalar-vector-tensor decomposition of the metric and scalar- 
vector decomposition of the lapse with 

d i hj = 0, d l hJJ = = d j hJT, d l nf = 0. (2.6) 
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We choose to define every derivative in the decomposition of the perturbation with a factor of 
a, such that the combination adx appears everywhere. 

Although apparently linear, the field perturbations contain in principle perturbations to 
all orders. For instance, 

tp = A^ 1 ) + AV (2) + 0(A 3 ) , A<1, (2.7) 

where A indicates the order in perturbations. Similar expansions hold for the other perturbations 
in (2.5). Note that the background value of iVj is zero (s, nf are of 0(A)), since we are describing 
perturbations on top of a FLRW spacetime. Moreover, having N(t) as a background for goo 
allows us to rescale time. For example, conformal time would be defined by setting N(t) — > 0(77), 
thus having N(t)dt — > a(rj)drj, with rj being conformal time. 

3 Gauge dependence and gauge invariant perturbations 

As mentioned above the gauge freedom in General Relativity corresponds to the fact that dif- 
ferent mappings can be chosen between the physical and background spacetime. Changing the 
mapping is then referred to as a gauge transformation. If x^ is the vector field associated with 
one mapping, then a different gauge choice x^ + transforms a quantity Q according to [7, 8] 

Q^e c zQ, (3.1) 

where is the Lie derivative along £. The 4-vector £ M contains all orders in perturbations, 
£ M = + 2^ 2 £(2) + 0(\ 3 ). £^ can be separated in a temporal and a spatial part, which can 
be written as 

e = (e,e) = U°,ld^+e {T A , (3.2) 

where is the transverse part of the spatial (di^ T ^ = 0). 

3.1 Gauge transformations of metric and field fluctuations 

We now study more precisely how the scalar and metric field in (2.1) transform under the 
gauge transformation (3.1). We consider the gauge transformations up to second order in A. 
For simplicity we only look at gauge transformations of the scalar degrees of freedom tp and 
Q, which transform under the temporal gauge parameter £°. We are not interested in second 
order perturbations of the other scalars n and V 2 s, because they are constraint fields that 
can be eliminated from the action by solving the constraint equations. Alternatively, they can 
be decoupled from scalar perturbations in a procedure [6] that constructs the gauge invariant 
action including gauge invariant constraints. Also, we neglect vector and tensor modes and 
spatial derivatives which are unimportant on long wavelengths. 

For the scalar field $ in (2.1) the gauge transformation (3.1) acts on the perturbation as 

<p _> ^ + 0A?£O + 1(^0 + + 2{p)R ^0 + ^ q.^ {3 3) 

Here the dotted derivative denotes a reparametrization invariant time derivative </> = jf^fi, such 
that it is easy to rescale time, for example conformal time by setting N(t) —> a(rj). Moreover, 
the spatial derivatives appear with a factor of a. At higher order we have only shown quadratic 
terms containing (temporal derivatives of) £°. Other terms include and/or spatial derivatives 
of (f or £°, which, as we mentioned, have been neglected. 
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Similarly, the metric tensor transforms under (3.1). If we consider the curvature perturba- 
tion £ then 



2( -)■ 2( + 2HNf + NC 



(3.4) 



Here H = a/ a. Note that £ only transforms under temporal gauge transformations t — > t + £° 
at linear order in perturbations. Also here we have not explicitly written other higher order 
terms that include and/or spatial derivatives of £ or £°. These terms appear when projecting 
out the gauge transformation of £. The precise transformations of £ and ip can be found, for 
example, in Refs. [7, 8] or [9] (the latter does not include the spatial gauge transformation £*). 

3.2 Constructing gauge invariant variables 

From Eqs. (3.3) and (3.4) a combination can be formed which is gauge invariant under temporal 
gauge transformations to first order. This combination is 



w 



H 

2^<p 



2C, 



(3.5) 



and is called the comoving curvature perturbation. This gauge invariant combination of scalar 
metric and field perturbations appears not to be unique, in the sense that it can be rescaled 
by any function depending on the background fields ((f), H). These background quantities are 
by construction fixed and do not induce additional gauge transformations. Thus there are in 
principle infinitely many gauge invariant combinations. However, a rescaling of a parameter by 
time dependent functions does not change any physical results. From that point of view the 
first order gauge invariant variable is unique. 

Note that we could have also made a gauge invariant combination of scalar metric or field 
perturbations with the perturbed constraint fields n or s which also transform under temporal 
gauge transformations. However, here we are only interested in dynamical gauge invariant 
variables. The constraint fields in the action (2.3) are nondynamical: they can be solved for and 
their solutions inserted back into the action. Thus, the only gauge invariant dynamical degree 
of freedom is the comoving curvature perturbation w. 

Although gauge invariant to first order in perturbations (governed by the small parameter 
A), to second order w changes under the gauge transformation 



= w + Ag (-w , 



H 



H4> H\ 



Hw + w 



m o +o(e,d^°,d iW ) 

(3.6) 



where the second order gauge transformation of w is A^ ^w, as in Ref. [10], with an additional 
subscript £ which indicates the gauge transformation involves terms oc £. w can be made gauge 
invariant to second order by adding quadratic perturbations to its definition. For example, we 
know that to first order £ changes under a gauge transformation as £ — > £ + HN^°. Therefore 
we can define 



W lr , 



w 



w 




1 



(C + wQ + —wC 



(3.7) 
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which is gauge invariant to second order 1 . This gauge invariant variable is related by a factor 
H/(f> to the gauge invariant field perturbation on uniform curvature hypersurfaces [8, 11], as 
one can construct it by fixing the vector field at each order such that £ = 0. Here, we have 
constructed the gauge invariant variable by demanding that the gauge transformation at each 
order is countered by appropriate terms. We emphasize that we do not want to use any gauge 
fixing in this work. Rather, we wish to maintain all fields and their perturbations, and eliminate 
any unphysical degrees of freedom by using gauge invariant variables. 

Alternatively we could have countered the temporal gauge dependence in Eq. (3.6) by 
quadratic terms in p which transform as tp — > p + <j)N^ . In that case it is useful to replace 
£ = -r<p — in the second order gauge transformation of w (3.6) such that 

\H4> ti J cf) 
= w + A^w . (3.8) 

A| piv here means the second order gauge transformation of w involving tp. Of course the gauge 

transformation here is the same as in Eq. (3.6), i.e. A| = A^w, because we have merely 
rewritten the gauge transformation in terms of p. Now we can define another variable 

Wf = w 

1 Hcj) 

= w — Ftplw, p] , (3.9) 

which is gauge invariant to second order as well (see footnote 1). In the literature [8, 11] this 
variable is constructed by fixing the gauge such that tp = at each order, and is therefore called 
the curvature perturbation on uniform field hypersurfaces. 

In principle we could have picked any combination of quadratic perturbations in £ and tp 
to balance the second order gauge dependence of w. The advantage of the above variables is 
that they reduce to the linear perturbations if one of the scalar perturbations is set to zero. 
For example, in the case where £ = 0, W v — > 2^<P, or when tp = 0, the other gauge invariant 
variable Wq —> — 2Q. As it turns out, this makes it very useful to find the gauge invariant action 
at third order in terms of these variables. We will discuss this in the next section. 

It is obvious that the two gauge invariant variables W v and are not equal at second 
order. Their difference can be expressed in terms of a gauge invariant second order part. For 
example, substituting the equality ( = Q<p — \w in the definition of (3.7), we find 

W ^ = W < + Jl[j I ) W C + ljj W C W C 

= W ( + Q(W C ,W ( ). (3.10) 

So, the difference between the two gauge invariant variables is quadratic in and its deriva- 
tives, and is therefore gauge invariant by itself. For clarification, at second order Wj? = W^, 

1 Here we have not written terms with vectors, tensors and spatial derivatives of w which are unimportant on 
superhorizon scales. The complete second order gauge invariant variables can be found in, for example, [8]. 
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so Eq. (3.10) gives a nonlinear relation between the different second order gauge invariant vari- 
ables. Both gauge invariant variables can be called 'physical' degrees of freedom, in the sense 
that they do not depend on the unphysical gauge degrees of freedom. The question is if the 
variables describe the same physics. One could imagine that, since the variables are not equal, 
the 2-point and 3-point functions are also different and this may give different results. Of course, 
in order to calculate the 2-point and 3-point functions and describe the dynamics of the gauge 
invariant variables, one should study the action for them, which is what we do next. 



4 The gauge invariant action for cosmological perturbations 
4.1 Gauge invariance at zeroth order 

The starting point is the Einstein frame action (2.1), which is manifestly covariant. Now one 
can insert the ADM metric (2.2) and one finds - up to boundary terms - the action (2.3). 
Although the general covariance is not manifest in this action, it is still present. We have merely 
decomposed the metric in separate parts, but as a whole it still transforms as a tensor. 

Now we insert perturbations on top of a fixed homogeneous, isotropic and expanding back- 
ground. If we consider the background alone, the action is 

= J d 3 xdtNa 3 |-3if 2 + ^ 2 - V(4>)\ . (4.1) 

This action is trivially covariant, in the sense that the background quantities transform at zeroth 
order under coordinate transformations, such that the background fields are fixed functions of 
the coordinates. In other words, if a = a(t) in one coordinate system, than a = a(t) in another 
coordinate system. 

The hamiltonian constraint, momentum equation and field equation for the background 
are found by varying the action with respect to the N, a and (j), respectively, 

3H 2 = + V{4>) 
2H = -<j) 2 

= + 3H(j) + V'((j)) . (4.2) 
It turns out that it is useful to define a variable z as 

« = |, («) 

such that the various slow-roll parameters can be written as, 

H 1 o 

e ^-JP = 2 Z 

Here the same definitions have been used as in [2]. These slow- roll parameters are very useful 
for finding the dominant contributions to n point functions from the action. 
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4.2 Gauge invariance of quadratic action 

The action to linear order in perturbations vanishes due to the classical background equations of 
motion. The first nontrivial action of perturbations is the second order action in perturbations. 
Usually one first eliminates the constraint fields N and iVj by solving for them and inserting 
their solutions back into the action. It is only necessary to do this to first order in perturbations 
since the second order solutions multiply the classical equations of motion. If we consider only 
the scalar fluctuations ( and <p, the resulting second order action looks schematically like 



where £ and ip are now linear in the perturbation parameter A, such that the action is of order 
A 2 . With the indication 0(( 2 ) we mean all terms of order £ 2 , including £ 2 , di(d l ( and 
possible mixings. A similar reasoning applies to the 0(ip 2 ) and O(Cv) terms. The explicit form 
(4.5) can be found in e.g. [12] or [6]. 

As was mentioned in Sec. 2 (and shown explicitly in Sec. 3), the perturbations £ and ip 
transform under linear gauge transformations, specifically under temporal gauge transforma- 
tions. Therefore, the variables £ and <p do not separately have a physical meaning, as they 
depend on our choice of coordinate system. The combination w in Eq. (3.5) is on the other 
hand gauge invariant, and if we express the schematic action (4.5) in terms of this w, one 
obtains [12] 



up to total derivative terms, z was defined in Eq. (4.3). Apart from w, the action also contains 
the transverse, traceless metric perturbation hy, in short the graviton, which is automatically 
gauge invariant to first order. 2 

The action (4.6) is now manifestly gauge invariant up to second order in perturbations: 
the variables are diffeomorphism invariant up to a first order in coordinate reparametrizations, 
and can therefore be considered physical, and the background fields are trivially gauge invariant. 
One can now do the usual steps of rescaling w to a variable v = jjW, the Mukhanov variable, 
such that the action becomes that of a harmonic oscillator with a time dependent mass. One can 
then quantize v with the usual expansion in creation and annihilation operators, solve for the 
mode functions in, for example, an inflationary background, and calculate the power spectrum 
for those perturbations that remain constant on superhorizon scales (which coincides with the 
comoving curvature perturbation w on superhorizon scales during slow-roll inflation). 

4.3 Gauge invariance of cubic action 

In order to study non-Gaussianities in the Cosmic Microwave Background radiation (CMB), n- 
point functions such as the primordial bispectrum or trispectrum must be derived. The 3-point 
functions are found from the third order action in cosmological perturbations. Of course, the 

2 Instead of solving for the constraint fields, it is also possible to keep the constraint fields and perturb them 
to second order [6]. After decoupling them from the dynamical degrees of freedom there are an additional 4 
non-dynamical degrees of freedom in the second order action, which are gauge invariant as well. Varying the 
action with respect to these degrees of freedom gives the first order solutions of the constraint fields. Thus, out 
of 11 degrees of freedom in the metric and scalar field, we are left with 3 dynamical and 4 constraint degrees of 
freedom. The 4 gauge degrees of freedom have thus been eliminated from the action. 




(4.5) 




(4.6) 
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physical bispectrum is only found when the perturbations in the (tree-level) cubic action are 
physical, i.e gauge invariant. 

Schematically, the third order action in £ and p takes the form 



£( 3 )(C, tp) = j d 3 xdtNa 3 {AO(( 3 ) + BO{Q 2 tp) + CO((p 2 ) + DO{p 3 )} 



(4.7) 



Again, one can find this action by dropping vector and tensor perturbations and solving the 
constraint equations to first order in perturbations. The third order solutions multiply the 
background equations of motion, and the second order solutions multiply the Hamiltonian and 
momentum constraint evaluated at first order. The AO{Q 3 ) and DO(ip 3 ) terms were first derived 
in Ref. [1], and for completeness we give the terms here explicitly. First the AO((, 3 ) terms up 
to temporal and spatial boundary terms: 



S®(Q=J d 3 xdtNa 3 -(C 2 + 2C-|) 



+ 



3C 



H' 

didjip didjijj 



(c4> 



a 



lz 2 
2H l 



a a 



(4i 



where 



Secondly the DO(p 3 ) terms: 



v 2 ( i v 



(4.9) 



5 , ( 3 )( (/J ) = / d 3 xdtNa 3 \ - 1 ztp~ip - -z I — - I tp - tp 



,2 . 1 . ( dup 



a a 



'didjx\ ( didjx 



+ 



-H 2 z 3 (3 — -z 2 ) — -zV" — -V'" 
8 v 2 ' 4 6 



1 



tp 6 + -z i Hp 2 Lp + -z z tp 



1 



.2,„2 



V 2 x 



(4.10) 



where 



(4.11) 



The AO((^ 3 ) terms have been derived for generalized scalar theories in Ref. [2], see also [3]. 
4.3.1 Manifest gauge invariance: cubic action for W v 

The action (4.7) does not appear to be covariant due to the gauge dependence of ( and ip. 
However, we know that the complete action is covariant, and should be covariant to this order 
as well. To make this more manifest, we can try to express the cubic action in terms of the 
linear gauge invariant variable w, since the third order terms in (4.7) only transform under linear 
gauge transformations. In doing so we have several possibilities. 

The first option is to eliminate all tp dependence in the action by replacing 



tp = ^zw + z(, 



(4.12) 
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where again z = <fr/H such that 




S( 3 ) = / d 3 xdtNa 3 \dO[ l-zw) ) + EO{( 3 ) + FO(C 2 w) + GO{Qw 2 ) \ = S§}{w)+S§l(C) • 



GI V^;-r^GDV 

(4.13) 

The first part S^j (w) now only depends on 0(w 3 ) terms, and is therefore explicitly gauge invari- 

ant at third order in A. The second part 5qq(C) is on the other hand gauge dependent, because 
it contains cubic terms £ 3 , C, 2 w and Qw 2 that cannot balance each other's gauge transformations. 
However, the original action (2.1) is diffeomorphism invariant, and this gauge dependence should 
somehow cancel. One should not forget that also the second order action changes under a gauge 
transformation. Although (4.6) is gauge invariant to first order in A, to second order in gauge 
transformations it is not. The gauge transformation (3.6) thus generates third order terms in 
A from the second order action (4.6). Following Rigopoulos [10], under a second order gauge 
transformation of w, 

S(2) H S (2) {w) + J d 3 xdt Na s ^ y^^ A| c w| , (4.14) 

where was defined in (3.6). In order to keep the general covariance at each order in 

perturbation theory, this second order gauge transformation of S^(w) must be balanced by 
appropriate linear gauge transformations of S^K This means that the gauge dependent terms 
in Eq. (4.13) must be proportional to the first order equations of motion, such that 

4d(0 = "/ d 3 xdtNa 3 |^^Mf c KC]| , (4-15) 

where F^[w, £] is defined in Eq. (3.7) 3 . Of course, the gauge dependent part 5*qd(C) is precisely 
such that it counters the second order gauge transformation of S^ 2 \w) in (4.14). Now it is 
straightforward to express the action in terms of the second order gauge invariant variable W^, 

S {2) + S® = I d 3 xdtNa 3 \z 2 ^-W 2 - \ (^) 2 | + S§(Wy • (4.16) 

Thus, we have found the third order action which is manifestly gauge invariant up to second 
order in gauge transformations. The dynamical scalar degree of freedom in this action is W v , 
where \zW^ is the field perturbation on uniform curvature hypersurfaces. We haven't mentioned 
explicitly what is S <y3 \W ip ), that is, what are the gauge invariant cubic vertices for needed 
for the calculations of non-Gaussianities. The easiest way to find this is to set Q = from the 
start. In that case W v coincides with ~tp and 

SgH W v)^ J d 3 xdtNa 3 {DO (if 3 )} . (4.17) 



3 This implies that the terms of C?(C 3 ) should drop out in the third order action (4.7) after the replacement of 
V (4.12). 
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Hence, S§}{W V ) can be found immediately after the replacement p> — > -zW^ in the DO(ip 3 ) 
terms of (4.7) [10]. Using these terms from [1] the result is 



-W*W V 



1 fdtW v 



1 



+ 16" 



V 2 



l H z* - -z 3 z 
8 2 



v 2 v 



1 „2^2 i 1 t7^5^ | 1 ij2^6^q 1„2^ „4 T /// ^ 3 T/ W 



+ + |#V(3 - ^) - -^V 



6 



w« aw; 



(4.18) 



This demonstrates the usefulness of the field perturbation on uniform curvature hypersurfaces 
\zW<p in combination with the uniform curvature gauge ( = 0. The cubic action (4.18) can 
be simplified further by doing partial integrations of the third and last lines and using the 
background equations of motion (4.2), 



S§}(W V ) = J d 3 xdt 
16 



No 3 1 



2 f v 2 



di_ 
V 2 



v 2 v 







z 




zH 


} 



(4.19) 



In Eq. (4.19) different orders in slow-roll can easily be distinguished. Using the slow- roll 
parameters (4.4) it is clear that the first line contains terms of order e 2 , whereas the second line 
is subleading in slow-roll. 

In the derivation of the gauge invariant action for W v we have set £ = 0. What happens 
if we would have taken <p = from the start? In that case the third order action (4.7) contains 
only terms AO(C 3 ), and Wtn is a nonlinear expression in £. Therefore the third order action must 
contain terms proportional to the first order equation of motion for £, which can be absorbed in 
the second order action by identifying W v (ip = 0). These terms proportional to the equation of 
motion were identified in [1]. It was shown that after the field redefinition £ — > ^W v (ip = 0), the 
cubic action for <p is obtained by replacing W v — > ^r'f - 111 our language of gauge invariance, we 
redefine Q to the gauge invariant variable W v (in the gauge ip = 0), then restore the dependence 
on (p in Wtn, and finally set Q = to get the action for 0(<p 3 ) terms. 

4.3.2 Manifest gauge invariance: cubic action for Wq 

In the previous section we found the manifestly gauge invariant cubic action in terms of W v , 
related to the field perturbation on uniform curvature hypersurfaces. Just as well we could 
have constructed the gauge invariant action for W^, the curvature perturbation on uniform field 
hypersurfaces. Starting point is again the schematic third order action (4.7) and we follow the 
same steps as in the previous part. Instead of eliminating ip in terms of w and £ as in Eq. (4.12), 
we eliminate £ 

^-cp- w . (4.20) 



As before, the third order action can then be separated in a gauge invariant part depending 
only on 0(w 3 ) terms, and a gauge non-invariant part depending on (p. Of course, this gauge 
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dependent part must again balance the gauge transformation of (4.14). Therefore 



S( 3 ) = J d 3 xdtNa 3 i^AO ^("^) j + E'0{^ 3 ) + F f O{ip 2 w) + G'(D(ipw 2 ) j 

= 4 3 I ) H + 5^(^) 

= sgiH " J d 3 xdtNa 3 | ^ F v [w, if] | . (4.21) 

Again we can do a redefinition of w to the second order gauge invariant Wt using (3.9) such 
that 

+ S ^ = I d 3 xdtNa 3 \z 2 \^-W 2 - \ (^) 2 | + S§(W ( ) • (4.22) 

Thus, following this procedure one obtains the manifestly gauge invariant action at third order 
expressed in terms of the gauge invariant variable Wq. The simplest way to find the gauge 

invariant vertices in Sqi (Wq) is now to set the field perturbation (p = 0. Then, the gauge 
invariant variable becomes W^{ip = 0) = — 2£, and 

S { q1(W c ) ^> J d 3 xdtNa 3 {AO (C 3 )} • (4.23) 

So, if we set (p = 0, the third order action for £ immediately gives the gauge invariant action 
in terms of the curvature perturbation on uniform field hypersurfaces after the replacement 
C — > — \Wq. Using the AO(Q 3 ) terms derived in [1, 2], 



where 



did did /V 2 -0 
a 2 a 2 V a 2 



+2 a**wtvvi (424) 

a a a z \ 



This demonstrates the convenience of working with the gauge invariant variable Wt in combi- 
nation with the gauge ip = 4 . 



5 Uniqueness of gauge invariant action 

In the previous section two third order actions for cosmological perturbations were derived which 
were manifestly gauge invariant up to second order in gauge transformations. The general trick 
is that the gauge dependent parts of the third order action could be absorbed in the second order 
action, which defined a gauge invariant variable. In Eq. (4.16) the gauge invariant cubic action 
was expressed in terms of W^, in Eq. (4.22) in terms of Wt. Comparing the gauge invariant 

4 In the derivation of (4.24) we have performed several partial integrations with respect to spatial derivatives. 
The corresponding boundary terms do not contribute to the bispectrum, contrary to the boundary terms for 
temporal partial integrations, which will be discussed next. 
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actions (4.16) and (4.22), we see that the parts quadratic in the gauge invariant variables are 
the same. Thus the tree-level propagator for W v is the same for Wq. The gauge invariant parts 
of the action which are cubic in the gauge invariant variables, Eqs. (4.19) and (4.24), appear 
not to be the same, 



S§(W„)^S^(W C ). 



(3), 



(5.1) 



This implies that the gauge invariant vertices for W v differ from those for Wq. On the other hand, 
the gauge invariant cubic vertices originate from the same action. This presents an opportunity 
to find out exactly how the gauge invariant actions differ. 

5.1 Non-linear transformations 

To illustrate what is the difference between gauge invariant actions for non-linearly related 
variables, let us take a general action for a second order gauge invariant variable Wx, 



S(W X ) 



d 3 xdtNa 3 -z 



1, 



-W 



x 



diWx\' [ , o(3 ), 



-) \+S K g{(Wx) 



(5.2) 



This variable is non-linearly related to another second order gauge invariant variable Wy, 

W X = Wy + Q(Wy, Wy) , (5.3) 

where the Q(Wy ,Wy) is a completely general function quadratic in Wy, which can include 
temporal and/or spatial derivatives of Wy. The action for Wy then becomes 



Y 



a 



+ / d 6 xdtNa A 



a 3 z 2 QW Y 



-ft 



Q-diWy 
a 



+ / d 3 xdtNa 3 { Q 



1 8S^ 2 \Wy) 



5Wy 



(5.4) 



where 8S^ 2 > /8Wy represents the equation of motion following from the quadratic action, i.e. 



1 5SW(W Y ) 
a 3 " 6Wy 



1 



a 3 z 2 W Y 



+ ^V 2 Wy 



These terms are related to the variation of the complete action as 

1 8S(W Y ) 1 8S^ 2 \Wy) 



8Wy 



8Wy 



+ 0(Wy) = 0. 



(5.5) 



(5.6) 



which vanishes by the variational principle. In the third order action the terms oc Q8S^ /8Wy 
therefore are zero at the cubic level upon inserting the solutions of the equations of motion (5.6), 
and consequently these terms do not contribute to the 3-point function at tree-level. However, 
contributions enter in the tree-level action at higher order that can contribute to the 2- and 
3-point function through quantum loops, which we do not consider here. 

Inspecting the tree-level cubic gauge invariant action for Wy (5.4) we see that the bulk 
action Sq{(Wy) coincides with the gauge invariant cubic action for Wx (5.2). In this sense 



- 13 - 



the bulk gauge invariant cubic action can be called unique. Thus, the evolution of the 3- 
point function is independent of the choice variables Wx or any nonlinearly related variable, 
characterized by Wy. 

Of course there are also boundary terms in the action for Wy (5.4). In the in-in or 
Schwinger-Keldysh formalism they can contribute to the 3-point function. In this formalism an 
expectation value may be defined as 



where p\W + {t- m ), W-{t- m )] is the density matrix at initial time tj n , which for a pure initial state 
equals \I'*[iy + (ti n )]^'[W A _(ti n )]. W± is here a second order gauge invariant variable on the + 
or — part of the complex in-in contour, which could be, for instance, Wx or Wy. In general, 
the operator 0{W{t)) depends on both W+ and W- fields. In the simple case of an equal-time 
3-point function, operator ordering is not important, one can drop the ± subscripts from W's 
and 0(W(t)) = W(xi,t)W(x2,t)W(x3,t). Coming back to the boundary terms, the spatial ones 
do not contribute to the 3-point function. On the other hand, the temporal boundary terms 
give in general a nonzero contribution. 

For cosmological correlation functions (5.7) the initial time is often taken at to = — oo. Strictly 
speaking, this is not the correct procedure. As t — > — oo the quantum field theory of gravity 
becomes strongly coupled (that is, in that limit the physical momenta k/a — > oo) and pertur- 
bation theory fails. In practice one can define an in-in expectation value by starting from some 
finite initial time to a t which perturbation theory is well defined. Thus boundary terms at to can 
contribute. One can define the initial state for, for instance, the gauge invariant variable Wx to 
be Gaussian. As a consequence the initial state for any other non-linearly related variable Wy 
is explicitly non-Gaussian. This initial non-Gaussianity is then evolved through the bulk action 

(3) 

Sqi[Wy). It is important to distinguish how much non-Gaussianity is dynamically generated 
from some Gaussian initial state, and how much comes from a potentially non-Gaussian initial 
state. For example, if one observes non-Gaussianity for the variable Wy, but one defines the 
initial state for Wx to be Gaussian, then some of the final non-Gaussianity finds its origin in 
the initial non-Gaussianity in Wy. 

Let us now discuss the contributions to the 3-point function coming from the temporal 
boundary terms at time t. There can be various types of boundary terms. Boundary terms of 
the type W 3 (t) cannot contribute to the 3-point function, because the 5- function in Eq. (5.7) 
forces the "+" and "— " vertices to be equal at time t. In fact, these boundary terms do not 
naturally appear after the transformation to a different gauge invariant variable, as can be 
seen in (5.4). Nonetheless these terms can appear after additional partial integrations of terms 
oc W 3 (t) and oc W 2 (t)W(t) in the bulk action. This was exactly done for the action of W v in 
going from Eq. (4.18) to (4.19), which justifies the use of (4.19) as the gauge invariant cubic 
action for W<p. 

One has to be careful with other boundary terms such as aW 2 W or f3WW 2 , because they 
give in general non-negligible contributions to the bispectrum. The reason is that these terms 
contain the canonical momentum Tlw, which has a nonvanishing commutation relation with W. 
In fact these type of boundary terms generate disconnected parts of the bispectrum. In the 




[VW + VW.]0(W(t))p[W + (t- m ) : W^t in )} 



x exp [i / dt' (L[W+{t')} - L[W-(t')]) 5 [W+(t) - W-(t)] 




(5.7) 
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general example above, it can be seen that the temporal boundary terms in the action for Wy 
are precisely of the form Sq(Wy) = f d 3 xHw Y {t)Q{t)- Consider now an example where 

Q{Wy, W y ) = a(t)W Y + P(t)W Y W Y . (5.8) 

Then using the expectation value as defined in Eq. (5.7) one can compute to lowest (tree- level) 
order that 5 

(Wx{x 1 )Wx{x 2 )Wx(x 3 )) = (Wy(x 1 )Wy(x 2 )W y (x 3 )} 

+ 2a(t) ((Wy(x 1 )Wy(x 2 ))(Wy(x 1 )Wy(x 3 )) + sym) 

+ m(iWY(xi)W Y (x 2 )){WY(x 1 )W Y (x 3 )) 

+ {Wy(x 1 )Wy(x 2 )){Wy(x 1 )W y (x 3 )) + sym) , (5.10) 

where sym stands for other cyclic contributions. In words: the 3-point function for the variable 
Wx is computed using the gauge invariant cubic vertices in Sqj(Wx), the result expressed on 
the left-hand side, but it is also directly related to the 3-point function for Wy, computed using 
gauge invariant vertices for W Y , plus additional disconnected parts coming from the boundary 
terms, which add up to the right-hand side of (5.10). Note that the relation between the 3-point 
functions in (5.10) can be immediately derived by inserting the non-linear relation between Wx 
and Wy (5.3) into the left-hand side of (5.10) and using Wick's theorem. Thus in order to 
compute the 3-point function of one gauge invariant variable, one can use the gauge invariant 
action for another variable (which may have a more convenient form) and add disconnected 
pieces according to the non-linear relation. 

Since time t in (5.10) is arbitrary, relation (5.10) holds also for t = t{ n , telling us how are 
initial non-Gaussianities in the variables Wx and Wy related. These variables also define some 
spacelike hypersurfaces an d £wy- For example, if the initial state on £wy is Gaussian, 

then the initial non-Gaussianity on Ti\\r Y will be given by the terms multiplying a and f3 in (5.10) 
evaluated at t = t- m . 

5.2 Practical example: different gauge invariant variables 

The general discussion in this section demonstrates that the cubic gauge invariant actions (5.2)- 
(5.4) for different 6 gauge invariant variables are related, in the sense that they both have the 
same, unique, bulk action, but they differ by boundary terms and terms proportional to the 
equation of motion. This is not always obvious. For example, the terms proportional to the 
equation of motion can be separated, partially integrated, and the remaining terms can be 
written such that it is not clear that they are total derivative terms. Therefore gauge invariant 

5 An alternative way to compute the 3-point correlator (5.10) is by making use of the so-called interaction 
picture [1, 13], where an expectation value to lowest order in perturbation theory is given by 

{0(W(t))) = -i f* dt'([0(W(t)),H int (t')]}. (5.9) 

J t 

If for the interaction Hamiltonian H[ n t one considers the part with the boundary terms in (5.4) only, which are 
of the form Hw(t)Q(t), then it is straightforward to find the disconnected pieces in (5.10) using the canonical 
commutation relations. 

6 Different here means non-linearly related. This in turn means that the gauge invariant variables coincide at 
the linear level, but differ at quadratic order. 
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actions for different, non-linearly related variables can appear very different. To illustrate this, 
let us now consider a practical example. The gauge invariant variables W v and Wt are non- 
linearly related as in (3.10), which is an example of (5.3). According to the above discussion, 
their 3-point functions should therefore be related 



{WJ Xl )WJx 2 )WJx 3 )) 



+ 



' {{W c {x l )W c {x 2 )}{W ( {x 1 )W c (x 3 ))+ S ym) + ... 



z 

2JH 



(5.11) 



where the terms of higher order in slow-roll have been neglected. Note that this relation can 
also be inverted to give the 3-point function of Wq in terms W v . This is purely based on the 
non-linear relation between the gauge invariant variables, but it should follow from the actions 
as well. The gauge invariant actions for W v and Wq coincide at the quadratic level, but the 
cubic actions Eq. (4.19) and (4.24) look very different at first sight. For instance, the action for 



whereas the action for Wt (4.24) seems to 



Wip (4.19) is of second order in slow-roll (e z 
be of zeroth order. Maldacena [1] showed that it is possible to relate the cubic action for (p in 
the C = g au g e with the cubic action for £ in the tp = gauge. The two actions differ by terms 
proportional to the linear equation of motion (which do not contribute to the tree-level action) 
and by some boundary terms. When translated to our language of gauge invariant variables, 
the cubic action for W v (defined on the C = hypersurface) can be related to Wq (defined on 
the (p = hypersurface), up to boundary terms. After many partial integrations of (4.24), the 
result is: 

2 



I 



,3 7 Na 3 (l 4 

d 3 xdt < -z 4 

8 [2 



1, 



1 



did a 



+ iT~ 



1 ss^Hwa ) 
+ Q(W W C )± Wl\ 



1 

2 

did* 



d l W c 



W ( + Wc 



V 2 



diW t 



-z 2 W c W 2 



zH 



5W, 



? (3) 

>gi 



(Wf) + j 



d 3 xdt 



C 

Na 3 



Q(W W ( ) 



1 5SW(W C ) 



8W C 



(5.12) 



This is precisely the form of the action predicted after insertion of the non-relation relation 
Wip = + Q(W^, Wf), as in (5.4), up to boundary terms. The gauge invariant vertices in the 

(3) 

bulk action Sq{ for coincide with those for W^. Thus, the cubic gauge invariant action for 
W v can be called the unique action which separates the different levels of slow-roll. A gauge 
invariant action in terms of any other gauge invariant variable, e.g. can be brought to this 
unique form after many partial integrations and extracting terms proportional to the equation of 
motion. The bispectrum for Wq can now be computed by making use of the partially integrated 
bulk action, including possible contributions coming from boundary terms which can give rise 
to disconnected contributions to 3-point functions, see Eq. (5.10)-(5.11). Alternatively, one can 
redefine the field to the non-linearly related W v and use the action for that gauge invariant 
variable. The 3-point function for Wq is then computed from Eq. (5.11). The non-linear relation 
between different gauge invariant variables prescribes what this field redefinition should be. 

5.3 Boundary terms, hypersurfaces and observations 

So far we have not shown the boundary terms in (5.12). They were explicitly computed for £ (or 
in <p = gauge) in Refs. [14] and [15] for scalar field Lagrangians which are a general function 
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<E> and its kinetic term. Both reach same conclusions: boundary terms with time derivatives of 
contribute to the bispectrum, and the dominant terms in the slow-roll approximation give 
exactly the same contribution as in Eq. (5.11), which is what one finds after a "field redefinition" 
to a non-linear variable in W^. 7 In our language this is nothing more than switching between 
different gauge invariant variables. 

One remark here is that the boundary terms in Refs. [14] and [15] do not disappear 
completely after redefining Wr to a new non-linear variable. On the other hand, Eq. (5.4) 
suggests that all boundary terms are incorporated after switching to a non- linear variable. This 
must be so because under a non-linear transformation 3-point functions of different variables 
are related as (5.10), irrespective of a specific action. When describing now different gauge 
invariant variables with their corresponding actions, the same relation of the bispectra should 
follow from the action, for both the bulk and the boundary. The origin of these additional 
boundary terms not removed by the field redefinition may reside in additional boundary terms 
of the quadratic action. The form of the quadratic action in (4.6) is only reached after several 
partial integrations, which generate additional second order boundary terms. Moreover, the 
original, unperturbed ADM action (2.3) also contains several spatial and temporal boundary 
terms that in principle contribute at every order. Together these boundary terms must add up 
to a gauge invariant second order boundary term, expressed in the linear w, since after all the 
original starting point is the covariant action (2.1). Now, under a non-linear field transformation 
these boundary terms will generate also cubic boundary terms, which may balance the extra 
boundary terms mentioned before. 

A different way to see this is to come back to the procedure of finding the gauge invariant 
cubic action, outlined in Sec. 4.3. Here the non-linear gauge transform of the second order bulk 
action was balanced by gauge dependent terms in the third order action, which are proportional 
to the linear equation of motion. This in turn defined a second order gauge invariant variable 
(Wtp or depending on the procedure). Similarly, also the second order boundary terms 
transform under non-linear gauge transformations. They can be written in a gauge invariant 
way after incorporating the gauge dependent third order boundary terms. 

In spite of these remarks the conclusion of [14] and [15] still stands: the dominant con- 
tribution to the bispectrum coming from the boundary terms for is taken into account by 
switching to a different gauge invariant variable Wm and using (5.11). This can be very useful. 
Ultimately we are interested in describing correlation functions of W{ (the second order gauge 
invariant comoving curvature perturbation) , because this is the field that is conserved on super- 
horizon scales and forms the initial fluctuations in the gravitational potential that are believed 
to form the seeds of structure formation and the observed temperature fluctuations in the CMB. 
The correct gauge invariant cubic action to use is therefore Eq. (4.24). However, this action 
does not clearly separate the dominant contributions in the slow-roll approximation. Before 
horizon crossing it is therefore much more useful to work with the non-linearly related variable 
Wm, for which it is straightforward to determine the dominant vertices. Note that one has to 
be careful concerning the initial state: if the initial state for is Gaussian at to, then for W v 
it will be automatically non-Gaussian. After horizon crossing the full action for should be 

7 In fact, in Ref. [15] the procedure is slightly different than stated above. After partial integrations Burrage 
et al. do not keep any boundary terms proportional W^Wq. The reason is precisely that these would generate 
disconnected parts of the 3-point function. The partial integrations performed are only the "allowed" ones: those 
that do not contribute to the bispectrum at all, or those that are slow-roll suppressed contributions. The field 
redefinition is slightly different in their work, and it does not coincide with the variable W<p (note: all computations 
are performed in the comoving, tp — 0, gauge.) 
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used, as this variable is conserved on superhorizon scales. 

Finally some remarks about non-Gaussianity and observations. As argued, the variable 
to use to calculate the primordial power spectrum and non-Gaussianity is W{. The CMB 
power spectrum, and possible non-Gaussianities there, are namely created by fluctuations in the 
gravitational potential. The universe reheats at a ip = (more generally 5p = 0) hypersurface 
where photons decouple. These photons carry the information of the gravitational fluctuations 
towards us. However, we observe these photons via satellites such as WMAP or Planck, that 
are to a good approximation freely- falling observers for which £ = 0, i.e. observations are made 
on zero curvature hypersurfaces. Choosing a certain gauge invariant variable to work with, is 
in essence nothing else but choosing the hypersurface. 8 We have seen that different gauge 
invariant variables are non-linearly related, and therefore their n-point functions are related via 
disconnected pieces. As a consequence, on one hypersurface perturbations may appear non- 
Gaussian, whereas they are in fact Gaussian on a different hypersurface, or vice versa. It is 
important to distinguish the amount of non-Gaussianity that is generated by the cosmological 
evolution, and non-Gaussianity that originates from the choice of hypersurface, either initially 
or at time of observation. 

6 Frame independent cosmological perturbations 

In the first part of this work we have discussed uniqueness of the gauge invariant action for 
cosmological perturbations with respect to the choice of gauge invariant variables. In this 
part we consider a different type of uniqueness for the action for cosmological perturbations, 
namely uniqueness with respect to different conformal frames. Examples are the Einstein frame, 
Jordan frame or string frame, which are related via field dependent conformal transformations 
of the metric and redefinitions of the scalar field. Our goal here is to show that at the level of 
perturbations the (cubic) action can be written in a unique form, which is independent of the 
Einstein or Jordan frame. 

6.1 Einstein frame and Jordan frame 

So far we have discussed perturbations of the Einstein frame action (2.1). It is also possible to 
study perturbations in the Jordan frame, where a non-minimal coupling between the scalar field 
and the Ricci scalar is present. The Jordan frame action is 

Sj = J d 4 x^gj\^-RjF{^j) - l -g^d^jd u <S>j - Vj{®j)\ . (6.1) 

All subscripts J indicate that the quantities are expressed in the Jordan frame. The function 
F(<&j) presents the general coupling between the Ricci scalar and the scalar field $j. Setting 
F($>j) = 1(= = l/(8nG)) takes us back to the minimally coupled case. An example of a 
model with nonminimal coupling is Higgs inflation [16, 17] where i ? ($j) = 1 + C^j, £ S> 1, and 
<£j is the Higgs field. 

8 Albeit the choice of a hypersurface seemingly relies on a gauge variant concept of setting e.g. ip = (comoving 
gauge), or £ = (zero curvature gauge), these surfaces have a well defined physical meaning in the following sense. 
Namely, choosing e.g. ip = completely fixes a gauge such that in that gauge —2C, can be associated with W(, and 
hence has a gauge invariant meaning. In that sense we can talk about a gauge invariant choice of hypersurfaces, 
Ew c and T*w v being just examples of two commonly used hypersurfaces. In this language, the question on which 
hypersurface one should perform calculations of non-Gaussianity becomes immaterial, as long as one knows on 
which hypersurface the observer measures, and how to relate the gauge invariant variables associated with the 
two hypersurfaces. 
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The Einstein frame and the Jordan frame are related via a combined conformal transfor- 
mation of the metric and a redefinition of the scalar field 

with 

n 2 = n 2 (*j) = f (6.3) 

Since these are just field redefinitions of the metric and scalar field, no physical information is 
expected to be lost in the frame transformation. This is what we refer to as physical equivalence 
of Jordan and Einstein frame. The physical equivalence is very useful, because it means we 
could obtain any results, such as the power spectrum, in the Jordan frame by transforming the 
well-known Einstein frame results using the above relations (6.2). Instead of dealing with the 
difficult nonminimal coupling, we merely have to deal with a modified potential. 

Although the physical equivalence between Jordan and Einstein frame, in the sense de- 
scribed above, is well established at the classical level, it is not obvious that it also holds at 
the level of quantum fluctuations. In order to solve this issue, one could derive the tree-level 
action for cosmological perturbations in the Jordan frame. As a first step, we again slice up our 
space-time using the ADM metric (2.2) and insert perturbations similar to Eqs. (2.5), 

9ij,J = aj(t) 2 {kj + hij,j(t,x)) 

$j = (f>j(t) + ipj(t,x) 

Nj = Nj(t)(l + nj(t,x)) 
N l j = aj{t)Nj{t)n i:J (t,x). (6.4) 

The background Jordan frame action (6.1) then becomes 

sf> = J d 3 xdtNja 3 j {-3H 2 F - 3HjF + i# - Vj^j)^ , (6.5) 

Here Hj = ^ = ^Mjdt± and ^ = d<j)j/(Njdt). 9 Moreover, F = F(<f>j), so F only depends 
on the background field <pj. The background equations in the Jordan frame are now obtained 
by varying the action with respect to Nj, aj and <j>j, 

3HjF = + Vj{4>j) - 3HjF 
2HjF = -<j) 2 j + HjF - F 

= h + ZHj4>j - 3(2^3 + Hj)F' + V'Mj) . (6.6) 

9 A dotted derivative on a Jordan frame quantity implies that it is a derivative with respect to Njdt, whereas 
a dotted derivative on an Einstein frame quantity implies it is with respect to Ndt. We will use this consistently 
throughout this work. 



^ 9liv,J 
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n 4 



Vj($j), 



(6.2) 
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Using Eqs. (6.2) we can find the explicit relations between background quantities in the Jordan 
and Einstein frame, 



N = 


MNj 


a = 




H = 






1 # h 


<P = 


$7 dtp j 




1 Q' 2 ■ 

& + (6.7) 

where Cl 2 = Q, 2 (tpj) = F(tpj). If we substitute these relations into the background action in the 
Einstein frame (2.1) or into the field equations (4.2), we immediately recover the background 
Jordan frame action (6.5) and its background field equations (6.6). This establishes the classical 
equivalence of Jordan and Einstein frame. 

6.2 Frame independence of second order action 

Next step is to find the second order action for cosmological perturbations in the Jordan frame 
by inserting (6.4) into (6.1). The nonminimal coupling between metric and scalar field is making 
this procedure much more complicated compared to the Einstein frame. It would be useful to 
establish also the frame equivalence at the level of perturbations, such that we can use the 
Einstein frame results and transform to the Jordan frame. 

Before we proceed, we notice that the metric and field perturbations are not the same in 
Jordan and Einstein frame. The conformal factor £l 2 = f2 2 (<3?j) can be perturbed to second 
order as 

0($j) = + fiVj + ^ / Vj, (6.8) 

where O' = d^l/dtpj. Using this and the general relations between Jordan and Einstein frame (6.2), 
we can write second order relations for ip and £ 

1 d 2 - 



2 



Obviously the Jordan and Einstein frame perturbations of the scalar field and the scalar part 
of the metric are not the same. Consequently, neither <pj nor £j are convenient variables to 
establish the quantum equivalence between the two frames, as these variables are inherently 
frame dependent. However, from the Einstein frame analysis we know that neither ip nor £ 
are physical perturbations, because they are gauge dependent. In the second order action for 
cosmological perturbations (4.6), we have seen that the only physical scalar perturbation is 
w, defined in (3.5). With the Jordan-Einstein frame relations for the background (6.7) and 
perturbations (6.9) at hand it is straightforward to show that at first order 

w = 2 ^<p-2( = 2^-<pj-2(j = wj. (6.10) 
9 <Pj 

The gauge invariant w in Einstein frame coincides with the gauge invariant wj in the Jordan 
frame to first order. Thus the physical comoving curvature perturbation is not only gauge 
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invariant, but frame independent as well. This frame independence to first order was first 
proven in Refs. [18, 19]. Using this one can immediately write down the second order Jordan 
frame action via a frame transformation from the second order Einstein frame action (4.6), which 
was done in, for instance, Ref. [20]. By making use of the relation 



z 



2 2 _ 1 ^ + 6^ 2 _ 1 , 

wwtTj. ; ti, 2 = w Zj ' (6 - n) 



{Hj + 



we can write 



Sf\wj) = J d\dtNj4\z^-^ - \ (^) 2 | • (6-12) 

Of course the second order action (6.12) is also obtained by starting with the original Jordan 
frame action (6.1) and doing the same steps as in the Einstein frame derivation. That is, use 
the ADM decomposition, solve for the constraint fields to first order, and write the action for 
physical perturbations. In Ref. [21] this was done without solving for the constraint fields, but by 
decoupling the constraint degrees of freedom from the dynamical degrees of freedom. The result 
for the scalar dynamical degree of freedom agrees with (6.12). This establishes the equivalence 
of frames at the level of quadratic perturbations. 

Note that when the action for perturbations is written in terms of gauge invariant and 
frame independent perturbations, there is no notion of separate frames. Although the Jordan 
frame and Einstein frame actions (6.1) and (2.1) look originally quite different (nonminimal or 
minimal coupling), at the level of perturbations the actions have a unique form. Therefore it 
does not make sense to talk about a preferred frame to describe perturbations. One can choose 
to start with an action for a certain nonminimally coupled scalar field, or for another minimally 
coupled scalar field (related to the first via the relation (6.2)) in a modified potential, and both 
give equivalent results. Of course, in the end one wants to express the results (such as n-point 
functions) in terms of variables in the defining frame. 

As an example, in Higgs inflation the defining frame is the Jordan frame, where the scalar 
field $j is identified with the Higgs field, the potential Vj is the Higgs potential and the field is 
coupled to the scalar curvature via F($>j) = 1 + C^j- One can directly compute perturbations 
in the Jordan frame, resulting in the manifestly gauge invariant quadratic action (6.12), but 
the derivation is difficult because of the nonminimal coupling. It is easier to transform to the 
Einstein frame first, and then use the standard result (4.6). Of course the scalar field in the 
Einstein frame is not the Higgs field, but is related to the Higgs field via a field redefinition, 
just like the Einstein frame metric. After finding the standard Einstein frame results one can 
transform back to the Jordan frame (the defining frame) in order to obtain (6.12). 

6.3 Frame independence of cubic action 

In the previous section it was shown that the comoving curvature perturbation in the Einstein 
frame coincides with the same perturbation in the Jordan frame at first order. At second order 
however w and wj are related as 



w = wj + 



'1 

(hj ' O, J d(j)/d(j)j \0 



Hj n'\ d 2 cj>/dct> 2 j _ n" fn' 



(6-13) 



Thus the first order gauge invariant variable w is not a convenient variable for comparing Jor- 
dan frame results with Einstein frame results at third order, since the variable itself is frame 



- 21 - 



dependent. Motivated by the results in the previous section, we could check if a second order 
gauge invariant variable is frame independent, in the sense that it has exactly the same form in 
the Jordan and Einstein frames. Indeed, it can be shown that the curvature perturbation on 
uniform field hypersurfaces, from Eq. (3.9), coincides with the same gauge invariant variable 
in the Jordan frame, 




2 1 ■ 

ip + — W(f 



Hj 2 1 . 

4>j 4>j 



(6.14) 



where we did not write the terms with vectors, tensors or spatial derivatives. The frame inde- 
pendence of was shown explicitly in [22], but other, more general, proofs exist as well [23-25]. 
This means that we can directly find the third order action in the Jordan frame from the third 
order Einstein frame action expressed in terms of W{, which can be used, for example, to find 
/nl for a nonminimally coupled theory [22, 23, 26]. We only have to replace the background 
Einstein frame quantities by corresponding Jordan frame quantities using (6.7). Thus from 
(4.24) we straightforwardly find 
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The 3-point function in for a non-minimally coupled theory is found in a similar way. Thus, one 
first transforms transforms to the Einstein frame where one can use previously computed results 
for the 3-point function for Wq, and this result can be expressed in Jordan frame quantities by 
going back to the Jordan frame. 

Of course we can also transform the partially integrated Einstein frame action for Wq (5.12), 
that shows the separation between different orders in slow-roll, to the Jordan frame. After the 
frame transformations we find 
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As explained in Sec. 5, the boundary terms are accounted for by performing a field redefinition 
to a new gauge invariant variable 



<p,j 



{zjWlj/ny 



(6.18) 



The 3-point function for W^j is then calculated from the 3-point function for W^j plus discon- 
nected parts in that variable 



+ 



Zj 



(W v ,j(xi)W v ,j(x2))<W ¥ ,,j(si)W v ,j(x 3 )) +sym + . . . , (6.19) 



where zj is defined in Eq. (6.11) and terms of higher order in slow- roll have been neglected. The 
variable W v j is of course the frame transformed W^, which is seen most easily after transforming 
both sides of the non-linear relation between gauge invariant variables (3.10). W^j is however 
not directly related to the field perturbation on uniform curvature hypersurfaces in the Jordan 
frame, since it does not reduce to in the gauge O = 0. In fact, the variable 



WJ 



(Cj + + -w(j 
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is precisely the gauge invariant variable in the Jordan frame that does this. The question is how 
to compute the 3-point function for the field perturbation on uniform curvature hypersurfaces 
in the Jordan frame: what is the action for W^ t j? It can be shown that Wu> j is non-linearly 
related to W v j as 
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(6.21) 



where terms higher order in slow-roll have been neglected. In the previous chapter 5 we have 
shown that the cubic actions for non-linearly related variables differ only by boundary terms, 
which can give disconnected contributions to the bispectrum. Thus the bispectrum for the field 
perturbation on uniform curvature hypersurfaces in the Jordan frame contains a connected part 
coming from the first lines of the action (6.17), plus a disconnected part from the non-linear 
relation (6.21). Alternatively, one could use the direct relation 



W v ,j = W C)J 




.22) 



by replacing all quantities by Jordan frame quantities in (3.10). Then one finds the connected 
part of the bispectrum for Wmj from (6.15), and disconnected pieces from (6.22). 

Now some words about the special situation when ipj = 0. In that case the conformal 
factor Q, 2 only has a background value and W{((p = 0) = Q = £■ Thus the cubic terms in (j not 
only directly provide the gauge invariant vertices, but can also be transformed to Einstein frame 
vertices, and vice versa. The third order Jordan frame action for £j was derived in [27], and 
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it was shown that it can be found from the Einstein frame action [1-3] in Ref. [28]. Moreover, 
one could imagine that at higher order in perturbations one can construct a gauge invariant 
variable which reduces to £j in the gauge <pj = 0, just as we did before for Wq to second order. 
When (pj = it now becomes almost trivial to show that the curvature perturbation is invariant 
under frame transformations. Thus, the curvature perturbation on uniform field hypersurfaces 
is frame independent to all orders [25]. 

Finally a remark about previous results found in Higgs inflation. Often computations are 
done in both frames, examples being quantum corrections of the Higgs potential [29-33], or the 
computations of the cut-off scale for which the theory becomes nonperturbative [34-39]. The 
Einstein frame results are then compared to direct Jordan frame computations by transforming 
them to the Jordan frame. More often than not, transformed Einstein frame results do not 
exactly agree with what is found in the Jordan frame. A recent example is a calculation of the 
field dependent cut-off in [39], which appears different in one frame or the other. However, the 
result of this section is that the Jordan frame action can be found directly from the Einstein frame 
via a field transformation. The most clear way to see this is that everything can be expressed in 
frame independent variables. Thus the cut-off scale should be the same whether you compute 
it directly in the Jordan frame, or via transformed Einstein frame results. The reason for the 
confusion and difference between results obtained in different frames in the references above is 
due to a non-covariant formalism, where the variables become frame dependent. For example, in 
Einstein frame computations often the nonperturbative cut-off scale is found from expanding the 
non-polynomial potential in powers of (p and neglecting metric fluctuations, ip is then a frame 
dependent variable. There is no confusion when using frame independent variables: quantum 
corrections or cut-off's computed directly in the Jordan frame are exactly the same when first 
computed in the Einstein frame and transformed back to the Jordan frame. The cut-off can 
therefore in principle be computed directly from (6.15), and from higher order generalisation of 
that action. 

7 Summary and conclusion 

In this work we have focused on scalar perturbations of a scalar field and the metric around an 
expanding background. These perturbations are generally gauge dependent. It is possible to 
construct gauge invariant variables by taking certain combinations of these perturbations. At 
second order in coordinate transformations there are in principle infinitely many gauge invari- 
ant variables, two specific examples being the comoving curvature perturbation and the field 
perturbation on uniform curvature hypersurfaces. These variables are related in a non- linear 
way. 

We have outlined the procedure for finding the gauge invariant actions for these variables 
at third order, and we have shown explicitly the gauge invariant cubic actions. Next we have 
demonstrated that, due to the non-linear relation between these variables, the cubic actions 
appear different, but actually their bulk actions are the same. They differ by boundary terms, 
which generically give disconnected contributions to the bispectrum. 

This brings us to the aspect of uniqueness. Once you pick a certain initial hypersurface, 
or equivalently choose a specific gauge invariant variable, for example the comoving curvature 
perturbation in a Gaussian state, then there is a unique action for this variable that evolves the 
initial (Gaussian) state and creates non-Gaussianity through the evolution. If the final hyper- 
surface is different then the initial, boundary terms must be taken into account which appear 
when switching to a different gauge invariant variable associated with the final hypersurface. 
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We also commented on initial non-Gaussianity: although for one gauge invariant variable the 
initial state is Gaussian, it is generically non-Gaussian for another variable. Some of the final 
non-Gaussianity for a certain variable therefore originates from some initial non-Gaussianity. 

Finally we discussed different conformally related frames. It was shown that the cubic 
action for gauge invariant perturbations in the Jordan frame can be obtained directly from the 
Einstein frame action. The trick is to identify the variable that has exactly the same form in 
either frame. The comoving curvature perturbation is such a frame independent cosmological 
perturbation. Thus the bispectrum for the comoving curvature perturbation can be found from 
the Einstein frame bispectrum by transforming the result to the Jordan frame. The bispectrum 
for another gauge invariant variable in the Jordan frame, such as the field perturbation on 
uniform curvature hypersurfaces, is then found from the non- linear relation between this variable 
and the comoving curvature perturbation. 

Frame independent cosmological perturbations can be a very useful tool in calculating 
quantum corrections for non-minimally coupled theories. An example is Higgs inflation, where 
often different results are found depending on which frame one uses for the computations due 
to non-covariance of the perturbed actions. These differences can be attributed to both frame 
dependence and gauge dependence of the formalism used. The issue is resolved when using 
frame independent, gauge invariant variables, since then the Jordan and Einstein frame actions 
have exactly the same form. 

The techniques used here can be quite straightforwardly generalised to compute higher 
order (quartic, etc.) gauge invariant tree level actions for scalar cosmological perturbations 
(needless to say, the choice of gauge invariant variables at cubic and higher orders is much richer) , 
and it is the subject of future study. Higher order gauge invariant actions are the necessary tool 
for gauge invariant calculations of quantum corrections to cosmological observables, the simplest 
one being the one loop correction to the equal time two point function for scalar cosmological 
perturbations. This work represents a step towards establishing a fully covariant framework for 
calculating quantum field theoretic correlators in curved space time backgrounds. Cosmology 
is a field where such correlators may in fact be observable. The problem of observables in 
cosmology is one of the most important unsolved problems in theoretical cosmology [40]. 
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